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■ demonstrate that in attractive channels the renormalization-group evolution of the 1+2 scattering 

. amplitude may be complicated by the presence of eigenvalues greater than unity in the kernel. We 

also show that these eigenvalues can be removed from the kernel by one subtraction, resulting in 
an equation which is renormalization-group invariant. A unique solution for 1+2 scattering phase 
Q ■ shifts is then obtained. We give an explicit demonstration of our procedure for both the case of 

^ , three spinless bosons and the case of the doublet channel in nd scattering. After the contribution 

of the two-body effective range is included in the effective field theory, it gives a good description 
;_j ■ of the nd doublet phase shifts below deuteron breakup threshold. 
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I. INTRODUCTION 



For over forty years now the three-nucleon problem has, with considerable success, been 
used as a testing ground for the nucleon-nucleon interaction. Usually the three-body equa- 
tions used are based on the Schrodinger equation or its implementation for scattering in the 
form of the Faddeev [l| equations with the two-body interaction being one of finite range. 
However, discussions of the case in which the range of the interaction is significantly less 
than the wavelengths of interest, i.e. kR <^ 1, also have a long history. There has been 
renewed interest in this case with the advent of Effective Field Theory (EFT) descriptions 
of few-nucleon systems at low energy j^, 3, 14 1 . In an EFT treatment of the problem of two- 
and three-body scattering at energies such that kR <^ 1 the two-body scattering problem 
requires renormalization since the leading-order two-body potential is a three-dimensional 
delta function. 

But in low-energy NN scattering k and R are not the only scales in the problem. The 
presence of a low-energy bound state in the A^A^ system — the deuteron — means that we 
must also account for the deuteron binding momentum 7 = a/— me^ ~ 1/a — with a ^ R 
the unnaturally-large NN scattering length — when we do an EFT analysis of this prob- 
lem. In technical terms the presence of an enhanced two-body scattering length — or equiv- 
alently a near-zero-energy bound state — means that there is a non-trivial fixed point in the 
renormalization-group evolution of this leading-order potential. As long as this is accounted 
for, and a power counting built around the scale hierarchy: 

~ 7 < 1/R, (1) 

a systematic EFT can be established and renormalized using a variety of regularization 



schemes , lid, lUl, ll^- A similar scale hiera rchy (and hence a similar EFT) 

governs the low-energy interactions of Helium-4 atoms It is also relevant to the 

physics of Bose-Einstein condensates, if the external rnagnetic field is adjusted such that the 
atoms (e.g. ^^Rb) are near a Feshbach resonance Q, Q]. 

As a first step in extending this EFT to heavier nuclei, the three-nucleon system was 
considered, and the Faddeev equations for the particular case of a zero-range interaction 
were solved. It was soon discovered that the leading-order (LO) EFT equation for the quar- 
tet (total angular momentum 3/2) channel yielded a unique solution [l^Q, while for the 
doublet (total angular momentum 1/2) channel the corresponding equation did not yield 



a unique solution — at least in the absence of three-body forces This could be 

simply understood on the grounds that in the quartet channel the effective interaction be- 
tween the neutron and the deuteron is repulsive as a result of the Pauli principle, and this 
ultimately means that the neutron and deuteron do not experience a zero-range interaction. 
In contrast, in the doublet channel the effective neutron-deuteron interaction is attractive 
and the full difficulties of the zero-range interaction manifest themselves. These difficul- 
ties were first elucidated by Thomas, who pointed out that — if two-body forces alone are 
employed — the nuclear force must have a finite range if the binding energy of nuclei is to be 
finite |2l|. 

The three-body scattering problem for zero-range interactions considered in the seminal 
wo.k of Beda,™ a„d coHabo.... Q U H fl Q wa. considered . .search 

that antedates Faddeev's landmark 1961 paper: by Skorniakov and Ter-Martirosian 22 1 
and by Danilov Q]. These authors found similar difficulties to Bedaque et al, and traced 
the non-uniqueness to the fact that in the asymptotic region this three-body equation for 
scattering reduces to a homogeneous equation whose solution can be added to the solution 
of the inhomogeneous equation with an arbitrary weighting — a point recently reiterated by 



Blankleider and Gegelia 
In their 1999 papers 



24 



13 



25 



26|. 



14| . Bedaque et al. introduced a three-body force into the 



leading-order three-body EFT equation, so as to obtain a unique solution for 1+2 phase 
shifts. They adjusted this force in order to reproduce the experimental 1+2 scattering 

14, The 



length. The energy dependence of the 1+2 phase shift was then predicted 
introduction of this three-body force is unexpected if naive dimensional analysis is used to 
estimate the size of various effects in the EFT, but it is apparently necessary if the equations 
are to yield sensible, unique predictions for physical observables. This also accords with the 
1995 paper of Adhikhari, Frederico, and Goldman, who pointed out that the divergences 
in the kernel of the Faddeev equations for a zero-range interaction may necessitate the 
introduction of a piece of three-body data so that these divergences can be renormalized 



away 



El- 



fBut see Refs. 



,120] for a conflicting view.) 
In an attempt to get some insight into alternative ways to establish a unique solution to 
the three-body scattering problem at leading order in the effective field theory, we try to 
bridge the gap between the Faddeev approach — in which the interaction has a finite r ang e — 
and the EFT formulation of this problem. In Sec. InJ we examine the Amado model |2^ for 
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the case of three spinless bosons. Here we look at scattering in which the interaction of an 
incident particle on a composite system of the other two is considered within the framework 



of the LaEfranffian for the Lee model 29 , 30|| . If three-body forces are neglected then the only 
difference between this approach and those at LO in the EFT of Refs. 

is that in the Lee model Lagrangian one may introduce a form factor that plays the role 
of a cut-off in the theory. In this way we can connect the LO EFT equations (without 
a three-body force) to those found in the Amado model, by taking the limit as the range 
of the interaction goes to zero. The resulting equation has a non-compact kernel unless 
a cutoff is imposed on the momentum integration. We then reproduce and reiterate the 



results of Refs. 



li 



23l |. demonstrating that the low-energy solution of the equation 



changes radically as the cutoff is varied. Using a renormalization-group analysis we trace 
this unreasonable cutoff dependence to the presence of eigenvalues equal to 1 in the kernel 
of the integral equation. 

In Section IIIII we use a subtraction originally developed by Hammer and Mehen 
to remove these eigenvalues. Our analysis of Sec. El then allows us to demonstrate that 
the subtracted three-body equation is renormalization-group- invariant. The subtraction of 
Ref. 0] was employed at a specific energy, and used experimental data from the three-body 
system to determine the half-off-shell behaviour of the 1-1-2- amplitude. Here we go further, 
and show that using low-energy two-body data plus just one piece of experimental data for 
the three-body system — the 1-1-2 scattering length — we can predict the low-energy three- 
body phase shifts. We do this by first solving the subtracted integral equation for the half-off- 
shell threshold l-|-2-amplitude. We then use this result to derive unique predictions for the 
full off-shell behaviour of the l-|-2-threshold-amplitude, and thence for the 1-1-2- amplitude 
at any energy. The equations derived in this way are equivalent to those of Bedaque et ai, 
but represent a reformulation of the problem in which only physical, renormalized quantities 
appear. In consequence the leading-order three-body force of Refs. HQ does not appear 
in our equations. Our single subtraction ultimately allows us to generate predictions for the 
energy-dependence of the 1-1-2 phase shifts at leading order in the EFT without the presence 
of an explicit three-body force. The subtraction does, though, require data from the three- 
body system (namely the l-|-2 scattering length) before other three-body observables can 
be predicted. 

In Sec. II VI we apply the formalism of Sees. HTVlHII to the — conceptually identical but tech- 
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nically more complicated — case of the doublet channel in nd scattering. Here, we compare 
the numerical solution to our once-subtracted equation with phase-shift data. In Sec. |3 
we consider higher-order corrections to the LO EFT and illustrate that the results from 
the EFT are in good agreement with the nd data below three-nucleon breakup threshold 
if the sub-leading (two-body) terms in the EFT expansion are adjusted so as to reproduce 
the asymptotic S-state normalization of deuterium. The resulting description of the dou- 
blet phase shifts is very good up to the deuteron breakup threshold. Finally in Sec. IVII we 
present some concluding remarks regarding the limitations of this method and discuss the 
convergence and usefulness of the EFT. 



II. THREE-BOSON SCATTERING AT LOW ENERGY 



Consider a system of three bosons at energies so low that the details of their interaction 
are not probed. Suppose, in addition, that two of the bosons can form a bound state — "the 
dimer", with binding energy — e^. In this section we will compute the amplitude for boson- 



dimer scattering in a low-energy effective theory. This problem has been studied for a. 



22l ]. and has recently been revisited in the context of EFT 



32 



most 



50 years 

In Sec. Ill Al we derive the Faddeev equations for this system, in order to establish the 
notation used elsewhere in the paper. We then show in Sec. IIIBI that this equation is 
not renormalization-group (RG) invariant, i.e. changing the regularization procedure used 
in the integral equation alters its physical predictions significantly. In particular, we will 
demonstrate that this lack of RG invariance is due to the existence of eigenvalues equal to 
one in the kernel of the integral equation. 



A. The Amado equations in the limit of zero-range interactions 

Consider a field theory of bosons A^, in which the two-boson bound state ( "dimer" ) D 
is included as an explicit degree of freedom. In this model the Lagrangian can be written 
as |24,Q: 

C = N^ lido + ^]N + D^AD-g[D^NN + DN^N^] . (2) 

n 

Here A is the bare inverse free propagator for the dimer. This is basically the Lee model [29[ 
for D ^ NN. Historically, in order to obtain a finite amplitude for boson-dimer scattering, 
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a regularization scheme has been invoked. This can be achieved either through the introduc- 
tion of a cut-off in all momentum integrals or by including a form factor in the interaction 
Lagrangian i.e. replacing g q{v), with p is the relative momentum of the two bosons 
m D ^ NN . The Amado model [2^ entails the second choice for the regularization. 
As a result the equation for boson-dimer scattering, after partial-wave expansion, takes the 
form 3 S4| 



X,(g, q'- E) = 2Z,iq, q'; E) + J dq" ^ 2^.(g, g", E) r [e^^^ - ^ j X,{q\ g'; E) . 



(3) 



with = E + irj, 7] a positive infinitesimal. The Born term q'; E) is the amplitude 
for one-boson exchange, and is given by 

+1 



A f 

Zi{q, q'; E) = - j dx 

-1 



9{K)g{Q) 



E-Mq^ + q'^ + q-q 



(4) 



where x = q ■ q', Pe is the Legendre function of order i and A = 1 for three identical bosons. 
In Eq. ^ the relative momenta of the pair in the vertices D ^ NN are given by 



K 



and Q 



(5) 



Note that the convention of Lovelace [36[ for the recoup 
by a factor of -1. In fact, as originally shown by Lovelace 



ing coefficient A differs from this 
Eqs. dni) and (jH) also govern 



nd scattering in the J = |; 5 = | channel, i.e. the quartet, but with A = —1/2 (in the 
convention used in this work). The details of the recoupling algebra for the bosonic, nd 
quartet, and nd = | channel, are given in Appendix 1X1 

The off-shell two-body A^A^ amplitude for this Lagrangian is of the form 



t{p,p';E)=g{p)T{E)g{p') 

where the dressed dimer propagator is given by 

S{E) 



(6) 



r{E) 



E-e, 

with ed the energy of the dimer bound-state (e^ < 0). The function S{E) is: 

dpp"^ g'^ip) 



(7) 



S{E) 



{E-Si){ea-^) 



(8) 
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and the residue of the A^A^ propagator at the dimer pole is then clearly S{ed)- 

At this stage the regularization function g{p) is present in the dressing of the dimer 
propagator t{E), and also in the one-nucleon exchange amplitude Zi{q,q'; E). In writing 
Eqs. (|7j) and (jH)), we have imposed the renormalization condition that the binding energy of 
the dimer takes its physical value. After this renormalization is performed the subtracted 
integral in Eq. (jH)) is finite and the cut-off function g{p) can be taken to be one. In that 
limit 



t{E) 



2 7 + y/-mE 



E 



TTm^ 



E-ed 



with 7 = ^J—med. The form of the two-body amplitude (jHl) is — in the limit g — > 1 — exactly 
that obtained at leading order in an effective field theory with short-range interactions 



alone 



11 



12 



This allows us to write the integral equation for 1-1-2 scattering in the limit g{p) 

S fE(+) - 



1 as 



q'- E) = 2Ze{q, q'; E) + / dg" g"^ 2Z,(g, g"; E) 
where, for 



X,(g",g';E), (10) 



0, the Born term 2Zo(g, g'; E) is given by 

g^ + g'^ + gg' — mE 



m 



2Zoiq,q';E)^Z{q,q'-E) = -\ — \n 

qq' 



g2 _|_ g/2 



qq' 



mE 



and 



\ Am , 



S{E;q) 



irm^ 



mE 



(11) 



(12) 



Here, we have included the factor of 2 resulting from the symmetrization for identical par- 
ticles in Eq. (jS)) in the definition of Z. These equations are identical to those employed 
in the EFT of Refs. 3, 3' although with different normalization for the amplitude Xi. 
The relationship between our conventions and those of Refs. 3, is elucidated further in 
Appendix |B| 

In the above analysis we considered the integral equation for the scattering amplitude 
or T-matrix. At any finite k the kernel of the integral equation has a pole, coming from 
the dimer propagator (see Eq. (fTn|) ). Furthermore, above the breakup threshold, the Born 
amplitude Z develops moving logarithmic singularities that need to be dealt with when this 
equation is solved numerically, for instance by performing a contour rotation that avoids 
these singularities. In this work we restrict the analysis to energies below the dimer-breakup 



7 



threshold, and so the logarithmic NNN cut is not an issue. In this energy domain a 
reformulation of Eq. (fTUI) that eliminates the D pole is useful. 
The reformulation involves writing the DN propagator: 

^ f - '^i {e-u- (13) 



We can then calculate the amplitude Xq of Eq. (fTUj) by first calculating the boson-dimer 
S-wave i^T-matrix using the integral equation: 

K{q, q'- E) = Z{q, q'; E)+V 1 dq" ^ Z{q, q"; E) Kiq", q'; E) . (14) 

E-ed-^ 

Below the three-body breakup threshold is a real symmetric matrix, and this equation is 
free of singularities. It is therefore numerically advantageous to solve this equation rather 
than Eq. ()in|l . and this is the approach we have used in generating our numerical results. 

The relationship of the phase shifts to the on-shell K-matrix K{k, k; E), is provided by 
first employing the relation: 

Xo(g, q'; E) = K{q, q'; E) - ^^K{q, k; E)S{E- k)Xo{k, q'; E). (15) 

Then, to determine the S-wave on-shell scattering amplitude and therefore the S-wave phase 
shift, we need to multiply the result from Eq. p5|) by the residue of the dimer propagator, 
i.e. define 

To(g, g'; E) = S'^^q, E) Xo(g, g'; E) S''\q' , E) . (16) 

The boson-dimer S-wave scattering phase shifts are then related to the amplitude Tq at the 
on-shell point by: 

where the on-shell momentum k is defined by the relation 

E-l^-U. (18) 
Using the relationships (fT^ and (fT7|) we find that: 

iT77 

K{k,k-E) = -—-iimS. (19) 
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Since we will be using the boson-dimer scattering length to renormalize our integral 
equation, we are particularly interested in the case = in Eq. (fTIH). Using the boson- 
dimer effective-range expansion: 

kcotS = -— + lr3p + ■ ■■ , (20) 
as 2 

with 03 and being respectively, the boson-dimer scattering length and effective range, we 
find: 

K{0,0-e,) = -^. (21) 
B. Renormalization-group invariance 

Now consider the convergence properties of the integrals in the equation ()14|1 . If q and 
q" are both large then Z behaves as 1/qq", while S scales as q". Therefore the convergence 
of the integral in ()14|) depends on the behaviour of X{q", k; E) at large q". Perturbation 
theory suggests that at large q" 

K{q\k-E)r^Z{q\k-E)r^l/q"\ (22) 

and thus the integral equation will be well-behaved without the need to impose any sort of 
regulator on the integral. However, this conclusion is erroneous. 

In fccLthe kernel ot d ha= inBnitely many eigenvalue, of order unity, as pointed out 
in Ref. [2^, and discussed in detail by Amado and Noble 37]. The argument of Amado and 
Noble may be heuristically paraphrased as follows. Calculating the trace of the kernel of 
Eq. (fT^ at E = ed we obtain: 

8A f dq" , ■ ' ' ^ 

tr kernel = — / m 

^ ' ?>Ti J q" 



U Jdq\ f3 q"^ + ^^ 

-,"2 I ^,2 



^ ^ ^ V 



3 , / 7 



I- '''' 



an integral which diverges, and does so logarithmically. Since the corrections proportional 
to (^-^^ do not affect the ultraviolet behaviour of this integral, if a cutoff A is imposed we 
have 

4 In 3 

tr(kernel) — ^ ln(A) (24) 
vSvr 

as A ^ CX3. Simple power-counting arguments demonstrate that the trace of all (positive- 
integer) powers of the kernel contain a logarithmic divergence too. Since it is also the case 



Qgarithmically with A — a 



that the largest eigenvalue of the kernel is finite the only way we can have 

tr(kernel°) ~ In (^^^ (25) 

for all n > 1 is if the number of eigenvalues larger than one grows 
result in accord with the analysis of Danilov [3| and Efimov 0, 

As we will discuss in more detail below, the presence of these order-one eigenvalues 
means that the perturbation theory estimate of K's large-momentum behaviour is incorrect. 
Consequently a cutoff must be imposed on the integral in (fT^ . as otherwise the integrals 
diverge. Once this is done the integral equation yields a unique solution for a fixed value of 
A. In this work we impose a sharp cutoff A, but choosing other, smoother cutoff functions 
does not alter the essence of the following argument. 

This regularization results in an equation: 

K^iq, q'; E) = Z{q, q'- E) + V dq" ^ Z{q, q"; E) ^ K^{q\ q'- E) , (26) 

^~^d- ^ 

where the presence of the cutoff has made K implicitly dependent on A. If the predictions 

of the theory are to be sensible we must have 

dKp,{q, q'-E) 

A ~ for q,q <t: A. (27) 

In other words, the low-energy predictions of our boson-dimer scattering calculation should 
not be affected by the imposition of a cutoff at a momentum scale far above the ones that 
are physically of interest. The demand (|77jl represents the renormalization group (RG) for 
this problem. 

Applying A^ to both sides of Eq. (j^Bj) we find that: 

A^{q,q';E) = A'Z{q,A-E) J^^''\, K^{A,q';E) 

+ V /rfg"g"^Z(g,g";E)^^^^A^(g",g';E). (28) 

For q,k <^ A 



4m 



Z{q,A;E)^j^; S{E;A)^A, (29) 

and, using the fact (which can be justified a postieri) that K{A, q'; i?) ~ ^ for A ^ q' the 
first term on the right hand side drops out of the integral equation ()28|) leaving: 

A^(g, k-E)=V Jdq" q"'Z{q, q"; E) ^ f ^^'f]^ ^^(^"' (30) 
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At this point it is tempting to argue that since Eq. has a trivial solution, it follows 
that K\ is RG invariant (up to terms of 0(1/A)). Such a conclusion is hasty, however. 
Non-trivial solutions to Eq. (|3Up exist if the kernel has an eigenvalue of one. In fact, the 
presence of eigenfunctions of this kernel corresponding to eigenvalue one has been proven 



rigorously in the asymptotic regime g ^ 7, by Danilov 



Here we repeat the analysis of Bedaque et al. 



which demonstrates the presence of 



23| and by Amado and Noble [S^. 



such an eigenvector. The impact of this "zero-mode" on the spectrum of the kernel has also 
been discussed by Gegelia and Blankleider 2^ Q and by Bedaque et al. 3^ . 

Consider E = e^, i.e. scattering at the threshold for the ND channel. Further, consider 
off-shell momenta q such that 7 ^ g ^ A. In this limit we may neglect terms proportional 
to 7 in S{E] q) and Z{q, q"; E), and so the integrand becomes scale invariant: 



dK\ 4A 1 



CXJ 

J dq" In 



q^ + q"^ + qq" 
q2 _|_ q//2 _ qq>f 



A 



dKA 

1a 



(g",0;e,). 



(31) 



Inspired by the absence of any physical scale under the integral we seek a power-law solution: 

A^(g,0;e,)~g-i (32) 



As shown in Refs. [13, 
equation: 



such a solution exists provided that s obeys the transcendental 
8A sin(f) 



1. 



(33) 



VSs cos (f ) 

If A = 1 this equation has complex roots s = ±isQ with sq = 1.0062. Thus in this case the 
RG equation has non-trivial solutions, which, provided A ^ g ^ 7, are of the form: 



f(q) = -Ccos ( soln ( - 



6 



(34) 



Note that Eq. (jHT|) also governs K{q, 0; e^), at least for A ^ q ^ 7, and so a similar analysis 
applied to K itself justifies the scaling K{q, 0; e^) ~ 1/q 



> q > 

HQ- 



For finite A we can examine the eigenvalue spectrum of the common kernel of Eqs. ()2f)|l 
and (jnni)- Doing this at = -v^ = and A = 1 for a variety of cutoffs yields the results 
shown in Fig. ^ As the cutoff increases there are more and more eigenvalues larger than 
one, with a new eigenvalue of one appearing each time the cutoff is increased by a factor of 
g7r/so — 22.7. This corresponds to an increasing number of bound states of the boson-dimer 
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system described by Eq. (l26|) . with the number of bound states growing as |37l |: 



I 

fi 



X2r 



"T 1 I I I I I 




10' 



1.0062, /A 
N = In - 

71 V 7^ 



(35) 



□ □ 



□ 



~~l — I — I I 1 1 1 1 II 1 — I — I I 1 1 1 1 II ^ 

o occaam — o ooooaBD — o ooooggami}) — o ooooodo — e-o 



1 — I I 1 1 1 1 1 




FIG. 1: The largest ten eigenvalues of the kernel of Eq. ^'AO^ . for the case A = 1, and a range of 
cutoffs. 

The presence of these eigenvalues which cross one as the cutoff is increased manifests itself 
as non-trivial cutoff dependence when Eq. (j26|l is solved. Some results found by solving this 
equation for the half-off-shell amplitude (again at zero energy) are displayed in Fig. |21 Here 
K{q, 0; td) ~ l/o', and so we have chosen to present results for the quantity: 

2nmS{E; q) 



a{q, k) 



-Xo{q,k;E). 



(36) 



This also aids comparison with the results of Ref. jUiJ, with which we are in complete 
agreement. The renormalization-group argument of this section ties the large changes in the 
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low-momentum amplitude seen in Fig. |21to the spectrum of the kernel of Eq. ()26p. via the 
concomitant strong-RG evolution of Ki^ at low momentum. 

3.0 I • • • • • • • 




FIG. 2: Results for a(g, 0), as defined in Eq. when various different cutoffs are imposed on the 
integral equation H14() . 

Note that in contrast to the case A = l, ifA = — | then the kernel of Eq. has no 
eigenvalues larger than one, and so the renormalization-group evolution of the K-matrix 
is smooth in that case. This means that as A ^ oo there is no cutoff-dependence in the 
p.edic«o,. fo. phase .h,f.-a. was see. .o be *e case „u,nerical,v Ref. Q where 
calculations of nd scattering in the quartet channel were performed in the EFT. (See also 

n 

Ref. J3j|.) The coupling at which the large-A RG evolution becomes non-trivial is the value 
at which Eq. (jHHjl first develops complex roots. This is [2^: 

(3T) 



13 



III. THE METHOD OF SUBTRACTION 



Thus, as it stands, if A > Ac, Eq. (fT^ needs additional renormalization before it can yield 



m 



14l | was to add a counterterm 



RG-invariant predictions. The solution proposed in Ref. 
to cancel the cutoff-dependence observed in Fig. |21 The three-body force introduced to 
renormalize the integral equation is not naively of the same order as the terms in the EFT 
Lagrangian but the analysis of Ref. 



m 



14l |. which has been recast in the previous section. 



shows that it is necessary for renormalization. The naive dimensional analysis estimate of 
the size of three-body forces is trumped by the presence of the shallow bound state in the 
two-body system, which is ultimately what leads to the Efimov spectrum shown in Fig. ^ 
Of course, as with any counterterm which removes cutoff-dependence in a quantum field 
theory, apiece of data is required to fix the value of the counterterm at a particular scale. 



a piece 

HQ 



In Ref. Iialllj the boson-dimer scattering 
More recently Blankleider and Gegelia 



eng th 03, was chosen for this purpose. 



21 



26[ have avoided introducing a three-body 



force in the leading-order three-body EFT equation by examining the solution of the ho- 
mogeneous equation and subtracting the oscillatory behavior. However, in their work no 
predictions for phase shifts were actually made. A subtraction technique for the three- 
body problem with zero-range forces was also suggested by Adhikari, Frederico, and Gold- 
man [27]. This technique was implemented in another integral equation with a non-compact 
kernel, that describing unregulated one-pion-exchange between two nucleons [4c|. However, 
the subtraction in Ref. ||4(| is performed at large negative energy, and involves demanding 
equivalence of the full and Born amplitudes at these energies. 



m 



In this work we suggest an approach which is equivalent to that used in Ref. 
but is formulated in an alternative fashion. Our procedure involves a subtraction of the 
on-shell amplitude at some — arbitrarily chosen but low — energy. The subtracted equation 
has a unique solution, which is, up to corrections suppressed by (p/A)^, RG-invariant. 



A. Subtraction at threshold 

Let us first consider the subtraction method applied to the integral equation for the half- 
off-shell 1 + 2 three boson threshold amplitude. We use information on the boson-dimer 
scattering length to fix the on-shell amplitude at the threshold. 
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Consider Eq. (fT^ for the half off-shell amplitude at E = e^, i.e. 

K{q, 0; erf) = Z{q, 0; ^d) - — j dq" Z{q, q"; e^) ^(e^; q") K{q", 0; e^) . (38) 





On the other hand, the on-shell amplitude at threshold should obey: 

A 



K{0, 0; e,) = ^ = Z{0, 0; e,) - ^ f dq" Z{0, q"; e,) S{e,; q") K{q", 0; e,) . (39) 



If we now subtract Eq. (jH^ from Eq. ()38p. we get an integral equation for the half-off-shell 
amplitude for which the input is the boson-dimer scattering length 03, in addition to the 
two-body data, which in lowest order is just the binding energy of the dimer. This equation 
is 

K{q, 0; e,) = ^ + AZ{q, 0; e,) - ^ f dq" AZ{q, q"; e,) S{q"; e,) K{q", 0; e^), (40) 
87 6 J 

' 

where 

AZ{q, q'; E) = Z{q, q'; E) - Z(0, g'; E) . (41) 



This equation (albeit in different notation) was first derived by Hammer and Mehen 31 1. 

In eq. ()4()|1 we have an integral equation in which the kernel goes to zero faster as q" 00 
than does that of the original integral equation. As a result we hope for a unique solution 
to Eq. pnjl . even if Eq. (fHHjl does not admit a unique solution. To establish this we need 
to prove that the amplitude 7^(5, 0; e^) is independent of the cut-off A, i.e. the solution is 
renormalization group invariant. Here we proceed as above, and differentiate the subtracted 
equation (jl0|) with respect to the cut-off A, to obtain 

^dK{^^ = -^AAZ(g,A;erf)S(A;e,)ir(A,0;erf) 

— ^ y ^^\<i^<i '^^d)S[q ;ed)A . (42) 



Once again, we consider 7 <^ g -C A, and in this regime: 

AZ(g,A;erf)~|^; S{e,-K)^K- (43) 

the inhomogeneous term in Eq. (jl^ goes to zero as q^ / h? for large A. Therefore, once again, 
in the limit K ^ 00 Eq. ()42|) is a homogeneous equation of the form 



dK{q,Q-ed) 4m ^ „ , dK{q",Q-ed) , . 

A = — J dq AZ{q,q ■,ed)S{q ■,ed)A ^ . (44) 
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It is now easy to show that the kernel of Eq. ()44p is negative definite. Thus, no matter 
how large we make A no eigenvalues of one can appear. This is demonstrated numerically 
in Fig. 121 where we plot the eigenvalues of the subtracted integral equation's kernel as a 
function of A/7. (Here we have chosen 037 = —2.) For the subtracted case, the kernel of the 
homogeneous equation, Eq. (jH)), has no eigenvalue close to one, thus there are no solutions 
to Eq. dO ^ - — Q • '• - ^ - X • • . ^nt of 
the cut-c 




9 r Q— QQCDQOQQcn rb 

50000 100000 150000 200000 

A/y 

FIG. 3: Eigenvalues of the kernel of the subtracted integral equation H14|) for various different 
cutoffs, for the case 037 = —2. 

In Fig. EJ we present, for several values of the cut-off A, the half-off-shell i^'-matrix at 
threshold 0; e^) that results from the subtracted equation. It is clear from the results 
that the solution of the subtracted equation is completely independent of A in the regime 
g <^ A, as anticipated from the RG argument above. In fact the cutoff can be numerically 
taken to infinity without any difficulty at all. 

Note that in the asymptotic regime 7 ^ g ^ A the subtracted equation (pn|l still has 
solutions 0; e^) of the form (jnH). These solutions ensure equality of the first piece of the 
integral in PU)) with the left-hand side of that equation, 0;erf). However, in contrast 
to the situation of the previous subsection, the solution of the subtracted equation in this 
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asymptotic regime is not scale invariant. It must still obey Eq. ()39|). since those pieces of 
Eq. (pn|) do not disappear when a solution of the form (jMj) is inserted. Thus — unlike the 
case of Eq. (fT^ — the asymptotic limit of (^01) is enough to determine the asymptotic phase 
6: 6 is fixed such that (jH^ is obeyed. 

3.5 f \ 




-9 -1 n 1 ? 4 R fi 

10 10 10 10 10 10 10 10 10 

q/y 

FIG. 4: Results for a{q,0) (see Eq. ()36() for definition) when various different cutoffs are imposed 
on the subtracted integral equation ((11]) . Here 037 was chosen to be -2. 

Thus our subtracted equation at threshold yields unique results for the half-off-shell 
amplitude without the need for an explicit three-body force. It also confirms that one piece 
of three-body experimental data is needed to properly renormalize the integral equation for 
the three-boson problem in the zero-range limit. 

B. The subtracted equation at any energy 

The above analysis was restricted to the amplitude at threshold and established that the 
solution of the subtracted equation is unique. The question now is: can we get the amplitude 
at any energy without any further subtractions? In other words: can we use the half-off-shell 
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amplitude at one energy and the original equation ()26p to obtain a RG-well-behaved K\ at 
all energies? 

To answer this, we need to write the on-shell amplitude at energy E in terms of the 
solution of the half off-shell amplitude at threshold. We do this in two stages. Rewriting 
Eq. (gOl) as 

A 

Avn r 

K{q, 0; e,) = K{0, 0; e^) + A[Z]{q, 0; ^d) - — J dq" A[Z]{q, q"; e,) S{q"; e,) K{q", 0; e,) , 



(45) 

and having determined the half-off-shell amplitude at threshold, we first need to determine 
the full-off-shell amplitude at threshold, i.e. K{q,q';ed). Before subtraction K{q,q';ed) 
satisfies the equation: 

A 

K{q, q'- e,) = Z{q, q'- ^d) - ^ j dq" Z{q, q"; e,) S{q"- e^) K{q", q'; e,) , (46) 



which has the original badly-behaved kernel of Eq. (j^Uj) . So, again we need to perform a 
subtractive renormalization. 

Since Z{q, q'; E) = Z{q', q; E), we have that K{0, q; e^) = K{q, 0; e^). But we know that 
K{0, q'; E), also satisfies the equation 

A 

477? f 

K{0, q'- e,) = Z(0, g'; ^d) - — j dq"Z{0, q"; e,) S{q"; e,) K{q", q- e,) . (47) 



We can now subtract this equation from the equation for the full off-shell amplitude — 
Eq. (gni)— to get: 








(48) 

This equation has the same kernel as Eq. ()45|) . and given that we have already determined 
K{0,q';ed) = K{q',0]ed), we can now determine the full off-shell amplitude at the elastic 
threshold. Numerical solution indeed confirms that the solution of Eq. ()48|1 is cutoff inde- 
pendent, and that the limit A — > oo can be taken. The resulting K{q, q'; ed) is also, by 
construction, real and symmetric. In this way we have established that the full, off-shell, 
amplitude at threshold can be determined with one subtraction, and therefore, given 03, we 
know the amplitude K{ed) = X{ed). 
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To derive the renormalized equation at any energy E for the amphtude X{E), we need 
to write the boson-dimer equation at the energy E, i.e. 



X{E) = Z{E) + Z{E) t{E) X{E) 



and the threshold equation 



X{ed)=Z{e,) + Z{e,)T{e,)X{e,) 



(49) 



(50) 



(Note that we will manipulate the equations for X, but the same manipulations can equally 
well be done with K.) These two equations can be written as: 

Z~\E) = X-\E) + r{E) (51) 
Z-\ed) = X-\td) + T{td) . (52) 

We now subtract Eq. from Eq. with the result that 

6[t] ^ t{E) - r(6,) = [Z-\E) - Z-\e,)\ - [x-\E) - X-\e,) 

= Z~\e,) [Z{e,) - Z{E)] Z-\E) - X-\e,) [X(e,) - X{E)] X'^E) . (53) 

Multiply this equation from the left by X(erf) and from the right by X{E), we get 

X{E) = X{e,) + [l + X{ed)r{e,)]6[Z]+X{ed)S[r]X{E) 

+ [1 + X(e,) r(erf)] 6[Z] t{E) X{E) , (54) 

where 

S[Z] = Z{E) - Z{ed) . (55) 

All integrals in the above equation have sufficient ultraviolet decay to be finite, with the 
possible exception of X{ed) T~{ed) 5[Z] t{E) X{E) which is a double integral. 

We now can write the above operator equation as an integral equation for the amplitude 
at a given energy E in terms of the fully-off-shell amplitude at threshold, X{q',q]ed), as 
input: 



X{q, k; E) = X{q, k; e^) + B{q, k; E) + J dq' q'^ Y{q, q'; E) X{q', k; E) 



(56) 



Where the second inhomogeneous term is 

A 



B{q,k-E) = 5[Z]{q,k-E) + J dq" q"' X{q,q";ed) r{ed;q") 6[Z]{q" ,k; E) , (57) 
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with 



S[Z]{q, k- E) = Z{q, k; E) - Z{q, k; e^) , (58) 
and t{E; q) = t{E — |^). Meanwhile the kernel of the integral equation is given by 

Yiq,q';E) = X{q,q';e,) S[T]iE; q') + S[Z]iq, q'; E) r(i?; g') 

A 

+ J dq"q"' X{q,q";e,) T{e,;q") 6[Z]{q" , q'; E) T{E;q') 



= X{q,q'-e,) 5[T]{E-q') + B{q,q'-E) T{E-q) , (59) 

with 

6[rmq) = T{E-q)-T{ea;q) . (60) 

In this way we can determine the half off-shell and from it, the on-shell amplitude, at any 
energy given the on-shell amplitude at one energy where the subtractive renormalization 
is done. Note that if 6Z = 0, i.e. the "potential" for the scattering equation is energy- 
independent, then 5 = and Y{E) = X{ed)6[T]{E). 

To test our procedure, we have calculated the boson-dimer phase shifts, for the case 
037 = 1.56. This value was chosen since models of the helium-4 dimer suggest a ratio of 
three-body and two-body scattering lengths of this size Figure El we plot kcot{5) 

against k (in units of 7) for five different cutoffs. Our results agree exactly with those 



reported in Ref. 
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14l |. In contrast to the figure presented in Ref. Q|, we see absolutely 
no cutoff dependence whatsoever in our results. No explicit three-body force is required to 
perform this renormalization. 



IV. NEUTRON-DEUTERON SCATTERING IN THE DOUBLET CHANNEL 

In the last section we formally developed the procedure for calculating the amplitude for 
1+2-scattering in the three-boson system at any energy, having renormalized the equation 
at threshold using the boson-dimer scattering length as input experimental data. The final 
results were identical for different cut-offs A. So far in our analysis we considered the three- 
boson problem in order to avoid the additional complication of coupled-channel integral 
equations. However, in order to establish the ability of the once-subtracted equations in 
EFT to reproduce experimental nd scattering data, we need to introduce the spin and 
isospin dependence of the nd scattering problem. The tensor interaction in the NN system 
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FIG. 5: Phase shifts for boson-dimer scattering in the case 037 = 1.56, for a number of different 
cutoffs. The curves are completely indistinguishable. 

is only manifest at higher order in the EFT without exphcit pions, and in this Section we will 
restrict our analysis to lowest order in this EFT, known as EFT(/), thus here we need only 
include the ^Sq and ^Si nucleon-nucleon channels. Since the nucleon-nucleon interaction in 
the ^Sq has anti-bound state, we can still write a dimer-like propagator in this channel, but 
now the subtraction point must be the energy of the anti-bound state. As a result we write 
the quasi-deuteron propagator as 



Ta[E) = — With a = s,t , 



(61) 



where s and t stand for the spin singlet and spin triplet nucleon-nucleon channels respectively, 
and: 

1 



mat 



= -- 



.ll- 

m ' 



= -7.88 MeV, 
7 = 45.71 MeV. 



(62) 
(63) 



The nd equations in the doublet channel are now a set of two coupled integral equations 
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in which the initial channel has the deuteron in the triplet (t), while intermediate states can 
have either a singlet (s) or triplet (t) NN pair with a spectator nucleon, all coupled to spin 
and iso-spin one half. These equations take the form 

Ku{q,q';E) = Zu{q,q'; E) +V J dq" q'" Zu{q,q"; E) n (e ~ ^] Ku{q",q';E) 

^ ^ 

^ / It2\ 

+V I dq"q"' ZUq,q";E)Ts [E - K ..{q" , q' ; E) (64) 

^ / //2 \ 

K,tiq,q';E) = Z,t{q, q'; E) + V J dq" q"^ Z^tiq, q"; E) n (e - ^\ Ku{q",q';E) 

^ ^ 

^ / //2 \ 

+V I dq" q'" Z,M. 4'\ E) Ts i^E - |-J K.^iq", q'; E) . (65) 



where v = |m is the reduced mass for the nd system, and the Born amplitude Z^p is given 
by 

m 

qq' [q^ + q'^ — qq' — ^E 
with the spin iso-spin factor matrix A given by: 

/ 



-7 ( ' T?\ \ 1 J + q''^ + qq' - ^e , 



(67) 



V 



1 -3 
-3 1 

It is not immediately apparent that the kernel of the coupled integral equations ()64|) and 
()65|) has the same problems as that of (j2SI)- By taking linear combinations of (jMj) and 
()65|) and looking in the asymptotic region we can perform an analysis akin to that used for 
Eq. (EHD Q- 

This shows that one subtraction is required to render the system 
well-behaved. Otherwise this kernel too, has eigenvalues which cross one as A is increased, 
and the RG-evolution of Ktt at low momenta will not be smooth. 

In this case the KuiS^, 0; e^), as given by Eq. (I7T|) . with 03 the doublet scattering length, 
is chosen for the subtraction. Experimentally j4lj : 

03 = 0.65 ± 0.04 fm. (68) 

We adopt the central value for 03. 

After the subtraction is performed the equations for the half-off-shell threshold amplitude 
become: 

X«(g,0;erf) = 7^*^(0, 0; e^) + AZii(g, 0; e^) 
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+ 1 dq"q"' AZuiq,q";E)TJe,-^ \ Ku{q",0;e,) 
^ ^ 

^ / //2 \ 

+ 1 AZi,(g,g";e,)r, - |^ 0; e,) (69) 

ir,t(g,0;ed) = Ku{0,0; e^) + [Zst{q,0; e^) - Ztt{0,0; e^)] 

^ / Il2\ 

+ 1 dq"q"' [Z,,(g,g";e,)-Z,,(0,g";e,)]rJe,-|- Kuiq",0;e,) 

^ / It2\ 

+ 1 dq"q"^ [Z,,(g,g";e,)-Zi,(0,g";e,)]rJed-|- ir,i(g", 0; e^) . 
^ ' 

(70) 

Once these equations have been solved for Ktt{q, 0; e^) and Kst{q, 0; e^) we can demand: 

Ku{q, 0; e,) = ^^^^(0, g; e,); 7^^,(0, g; e,) = K^tiq, 0; e,), (71) 

and so arrive at two sets of two coupled equations apiece. These four equations determine 
the fuUy-off-shell threshold nd scattering amplitude. The first pair is: 

Ktt{q,q']ed) = Ku{0,q';ed) + AZtt{q,q';ed) 

+ J dq" g"^ AZttiq, q\ td) '^\^'^~\~ \ ^ttiq", q\ e^) 
^ ' 

+ I rfg" AZt,(g, g"; e,) r, L - ^"j iT.^ (g", g'; e,) (72) 
^ ' 

Kst{q,q';ed) = Ktt{0,q';ed) + [Zst{q,q';ed) - Ztt{0,q';ed)] 

^ / Il2\ 

+ 1 dq"q"' [Z,,(g,g";e,)-Z,,(0,g";e,)]rje,-|- Kuiq" , q' ; ed) 
^ ^ 

^ / It2\ 

+ 1 dq"q"^ [Z,,(g,g";ed)-Zi,(0,g";e,)]rJerf-|- ir,i(g", g'; e^) , 

(73) 

which have exactly the same kernel as ()69|) and ()7U|) . but different driving terms. 

The second set of subtracted equations describes the (unphysical) amplitudes Kts and Kgs 
at threshold. The unsubtracted versions of these equations are given by a simple extension 
of Eqs. (jMjl and After subtraction the equations are: 

Ktsiq,q']ed) = Kts{0,q']ed) + AZts{q,q']ed) 
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^ / Il2\ 

+ j dq"q"^AZu{q,q";e,)Tt U - K^q" , q'; e^) 
^ ^ 

^ / //2 \ 

+ I dq" q"'AZUq, q"; ed)r, l^e, - |- j K^q", q'; e,) (74) 



Kss{q,q';ed) = Kts{0,q';ed) + Zss{q,q';ed) - Zts{0,q';ed) 



+ j dq" q"'[Z4q, q"; e,) - Zu{0, q"; e,)]n L - ^) K,,{q", q'; e,) 

A / //2 \ 

+ I dg" g"; e,) - ^^.(g, q"; e,)]r, (^e, - |- j K^q", <1 



ed). 
(75) 

Note that imposing ()7ip to perform the subtraction on the set of four original integral 
equations (written in matrix form in Appendix El leads to a symmetric result for the 2 x 
2 matrix form of the threshold amplitude. 

Now we write the original, unsubtracted, equations in operator form, as: 

K(E) = Z(E) + Z(E)r(E)K(E), (76) 

with K the 2x2 matrix: 

K=|^"^"|; (77) 



n 



_//2 ^ 



and Z the 2x2 matrix defined by Eq. 

We can then perform the formal manipulations that lead to Eqs. ()56|) - ()60|) . except that 
now all quantities are 2x2 matrices in channel space, and thus the final integral equation 
to be solved is, in matrix form, but with the momentum-dependence made explicit: 

A 

Xiq,k-E) = X{q,k-ed) + B{q,k;E) + J dq" q'" Y{q, q"; E) X{q" , k; E), (79) 



with: 



B(E) = 6[Z] + Xied)Tied)6[Z] (80) 
Y{E) = X{ed)6[T]+B{E)r{E), (81) 
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where the meaning of the energy-difference operator 5 is exactly as in the boson case of the 
previous section. 

Applying these equations to scattering in the nd doublet channel below the nnp breakup 

threshold yields the phase shifts shown in Fig. El At almost all energies shown these agree 

with the leading-order results published in Ref. [32| at the 1% level. Once again, there is no 

cutoff dependence, once the doublet scattering length is used to subtractively renormalize 
the eq " ■ - - - , „ , , . „ . [~| 
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FIG. 6: Phase shifts for neutron-deuteron scattering at low energies, at leading order in the nucleon- 
nucleon EFT without explicit pions, for a variety of different cutoffs. The curves for different cutoffs 
are indistinguishable. The triangle is the central value of the scattering length measurement of 
Ref. ^41^. The circles are the results of the van Oers-Seagrave phase-shift analysis j^l, and the 
squares represent a recent calculation of these phase shifts by Kievsky et al. [4^. 
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V. THE ND DOUBLET CHANNEL BEYOND LEADING ORDER 



In this section we discuss calculations of nd doublet scattering which go beyond the 



leading-order calculation of the previous section. That computation emp, 
grangian extended to the NN system. This Lagrangian is equivalent to 



oyed the La- 



£(0) = c = NUido + 



2M 



(82) 



with the spin-isopsin projector which restricts the interactions to the ^Si or ^Sq channel, 
as appropriate: 

P: = T-2r,; = ^a^a, r^, (83) 

and ^Co"2i is the leading-order contact interaction in these channels. The subscript o,-i on 
this coefficient indicates that it appears in front of an interaction which has no derivatives, 
but that it scales as Q^^, with the enhancement over its naive dimensional analsysis scal- 
ing being due to t he p resence of the unnaturally-large scattering lengths in the two-body 



system 



due to tne_pr 



The higher-order calculation we report on in this section requires the insertion of higher- 



derivative four-nucleon operators. The analysis of Refs 



n 



lOl Il2j indicates that the first 
additional piece of the EFT Lagrangian which must be considered is ^ 



£(1) ^ _ ^ ^cg(Ar^P°Ar)t(Ar^P°Ar) - [{N^P^NyiN^Ol'^N) +h.c. 



(84) 



a=s,t 



and the Hermitian, two- derivative three-component, operator O is defined by: 



O: 



2,o 



(85) 



Here the effect of the two-derivative operators on the A^A^ amplitude is suppressed by one 
power of the small parameter jR [R/ag in the ^Sq case) relative to the leading-order EFT 
amplitude. Also appearing in C^^^ is a small correction to Cq, denoted by ^Co,o: "small" 
because ^Co,o is down by 7P relative to ^Co_i 0,0, 1^ 1^ - 

Thus the effects of the terms in C^^'' on the A^A^ amplitude can be calculated in per- 
turbation theory. '^€2 can then be chosen so as to reproduce the asymptotic S-state 



normalization of deuterium 



4J, and ^C^'q'^ adjusted in such a way that double-pole term 



which would otherwise appear in the A^A^ EFT amplitude is removed. This produces a 
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next-to-leading order (NLO) ^Si A^A^ amplitude P|: 

Tip) = — * 
vrm 



'J + ip 27 



(86) 



with p = y/ ME, and where is the residue of the '^Si T-matrix at the deuteron pole p = ij. 
This amplitude is easily seen to be a re-expanded version of the effective-range-theory ^Si 
amplitude: 

^(P) = 1 f L 2^ ^ ■ ^ (87) 

vrm 7 - ^pt{p^ + 7^) + tp 

where pt is the A^A^ ^Si effective range, which is of order the range of the A^A^ interaction: 

Pt ~ R. The re-expansion is thus in the small parameters pt'y and ptp, but with 7 treated 

as being of the same size as p. By making such an identification we determine that: 

1 



' 1 

1 -7Pi 

Zt is also related to the asymptotic S-state normalization of deuterium, As |4^, : 



(88) 



Al = 2-fZt. (89) 
Using the Nijmegen PSA value for As, As = 0.8845 fm"^/^ we obtain: 

Zt = 1.686, (90) 



which agrees with the result obtained from Eq. (j88j) to three significant figures. 

To summarize, the coefficients in the NLO EFT Lagrangian may be chosen such that 
the amplitude in the '^Si channel has a deuteron pole with the experimental binding energy 
and the "experimental" asymptotic S-state normalization. Also present in the NLO A^A^ 
^Si amplitude is a constant piece, which is proportional to pt- Here we wish only to assess 
the impact of higher-order terms on the nd phase shifts, and thus, we will perform a partial 
NLO calculation of the nd phase shifts below breakup threshold, dropping the non-pole term 
in Eq. (18511 . Work on complete higher-order calculations within our subtractive framework 
is in progress, and these numerical studies, as well as prior results by other authors 1^ 
indicate that including the constant term of Eq. (j85|) has little effect on nd phase shifts 
below nd breakup threshold js^. 

Similar results follow for the NLO ^Sq amplitude, and there: 

Zs = . \ , (91) 
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r, = 2.73 fm |46| being the effective range in this channel. This results in much smaller 
NLO corrections from this channel, since Tg/as is only of order 10%. 

Thus, to perform our (partial) next-to-leading-order calculation for nd scattering the only 
changes to the amplitude which are necessary are the multiplication of and by factors 
Zt and Zs- The subtractive procedure developed above is not affected by the inclusion of 
these factors: the only changes necessary in the above equations are the replacements 



St StZt, 



SgZg. 



(92) 



Making these replacements we obtain the results shown in Fig. [7| Once again the result 
is cutoff independent. It a gree s remarkably well with the sophisticated potential- model 



calculation of Kievsky et al. 
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FIG. 7: Phase shifts for neutron-deuteron scattering at low energies, at leading order (dash-dotted) 
and next-to-leading order (solid) in the nucleon-nucleon EFT without explicit pions. The curves 
for different cutoffs are indistinguishable. Data are as described in Fig. El 



These results are a (partial) NLO calculation of the nd phase shifts below breakup. They 
differ from those of Refs. |32, |33| since in those works the authors chose to adjust the NLO 
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coefficients in the EFT Lagrangian to reproduce pt exactly, and so only obtained Zt (or 
equivalently As) approximately. The difference between this "p-parameterization" and our 
"Z-parameterization" is a higher-order effect, and the magnitude of the discrepancy between 
the results of Fig. [7| and those of Refs. 3^ |33| is consistent with an effect of order {pt'lf. 
i.e. two orders beyond leading. Comparison of our numerical results with those of Ref. [33| 
indicates that if we adopt the p-parameterization the agreement is better than 1% js^. 

VI. CONCLUSIONS 

The integral equation which describes 1+2-scattering in the effective field theory with 
short-range interactions alone does not yield an RG-invariant low-energy amplitude. By 
performing an RG-analysis of the integral equation for this process we traced this poor RG 
behaviour to the presence of eigenvalues of order one in the kernel of the integral equation. 
One subtraction removes these eigenvalues from the spectrum of the kernel and renders it 
negative definite (at threshold). Imposing Hermiticity and employing a series of resolvent 
identities we can use this single subtraction at the 1+2-threshold to generate predictions for 
phase shifts at finite energies. Note that although here we have only computed phase shifts 
below l-|-2-breakup threshold our subtractive technique is easily extended to include energies 
above the three-body threshold. The only complication is the technical one of dealing with 
the logarithmic branch cuts that appear in the kernel of the integral equation at this energy. 



The equations we have developed are equivalent to the equations of Bedaque et al. jl3|, , 
and may be obtained from those equations by algebraic manipulations. The distinguishing 
feature of our formulation is that the equations are subtractively renormalized, i.e. only 
physical quantities appear in them, and any regulator can be employed. This would appear 
to make this formulation especially useful for higher-order computations in the nd system. It 
also provides particular emphasis to the point that - as in the case with all bare parameters in 
field-theoretic Lagrangian - the three-body force which appears in the equations of Bedaque 
et al. is not an observable. 

Thus one piece of three-body experimental data is needed in order to renormalize the 
three-body equations for zero-range forces. For this piece of data we choose the 1+2 scatter- 
ing length. Its value can be incorporated into the EFT description of the three-body system 
either via a counterterm, as in Refs. HQQ or, as done here, by a subtraction of the 
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badly-behaved integral equation. The renormalization of the equation after the inclusion of 
this single piece of three-body data provides a simple, model-independent, explanation for 
well-known features of the three-nucleon system such as the Phillips line 
tates the systematization of predictions made by Efimov for such systems 

Finally, we performed a partial treatment of next-to-leading order corrections to the dou- 
blet nd phase shifts in the EFT. We found that adopting coefficients in the NLO EFT(/) La- 
grangian that give the correct deuteron binding energy and asymptotic S-state normalization 
results in excellent reproduction of potential-model 5 = | nd phase shifts below deuteron 
breakup threshold. Our results suggest that — to a very good level of approximation — these 
nd phase shifts are determined by four numbers from the two-body system, e^. As, and the 
^So scattering length and effective range, together with the crucial one piece of data from 
the three-body system: the nd doublet scattering length. 



Acknowledgments 

D. R. P. thanks Silas Beane, Paulo Bedaque, Harald GrieBhammer, and Hans- Werner 
Hammer for useful discussions. We also thank Hans- Werner Hammer and Harald GrieBham- 
mer for providing us with their results, and with phase-shift data for nd scattering in the 
doublet channel. D. R. P. is grateful for the hospitality of Flinders University, where much 
of this work was done, and that of the Institute for Nuclear Theory, where it was com- 
pleted. The work of D. R. P. is supported by the United States Department of Energy 
under DE-FG02-93ER40756 and DE-FG02-02ER41218. The work of I. R. A. is supported 
by the Australian Research Council. 



APPENDIX A: SPIN-ISOSPIN FACTORS FOR THE AMADO EQUATION 

In this appendix we derive the spin-isospin factors for the Amado model for: (i) three 
bosons; (ii) nd quartet; (iii) nd doublet. In the latter two cases we restrict our analysis to 
S-waves only. 

The Amado equation can be written in operator form as 

X = 2Z + 2ZT X 



30 



where 

= (1 - 5a(3){ga\Go{E)\gp) 

with Gq{E) = {E — Ho)^^. Note: This differs from Lovelace ^] by (—1) due to a different 
definition of Go- 

This Zai3 can be written after partial wave expansion as 

where Xap is the product of a spin factor A^^ and an isospin factor A^^, i.e. 

Kfs = Af^ A^^ . 



with 

Set "^7 Sfj 
■5/3 "S* 5*0, 

where Sa, Sj3 and are the spin of the three particles, and Sa is the total spin of the pair 
(P'j). This expression can also be used to calculate the isospin factor A^^. 

For three boson the spins and the isospin of all three particles is zero. In this case we 
have only one channel, i.e. Xap = A, and therefore the spin-isospin factor is one, i.e. 

Xap = 1, three bosons . 

For nd scattering all spins and isospins are |. For the quartet state we have only one 
channel with 5' = | and I = \- The spin and isospin of the pair is 1 and (the quantum 
numbers of the deuteron) respectively. In this case 

A^ - 1 A^ - -- 



and therefore 

1 
2 



Aq,/3 = — - , nd quartet . 



Finally for the case of nd doublet we have two channels. They correspond to the pair of 
nucleons being in either t = {Sa = l,ta = 0) the deuteron, or in s = {Sa = 0,ta = 1) the 
singlet. In this case S = I = ^. The spin isospin factors are 

A^- ^ A^- ^ ^ A - ^ 
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and 



or 



AS _ _^ = ^ 

~ 2 ' 2 



3 

1 



A 



a/3 



1 -3 
-3 1 



nd doublet 



APPENDIX B: CONNECTION TO EQUATIONS OF BEDAQUE ET AL. 

The equation of Refs. 



3, 



is, in the case of no three-body force: 



2A f 

a{k,p) = XM{k,p; fc) H / dqM{q,p; k) 

vr J 



qi _ _ 



a{k,q), 



with 



and: 



3 \ V 4 pq \q'^ — qp + p-^ — mh 



X 



1 for three bosons, 
— I for the nd quartet channel. 
Here the relation to the phase shifts is given simply by: 

1 ' 



Re 



k cot S. 



^a{k, k) ^ 

To make the connection to Eq. pO|) first observe that: 



S{E;q) = -^i^ + J^q^-mE\, 



4- 



and then define X{p,k : E) such that: 

4 



we then find: 



3 I ^ + \/~4~ I ^(P'^S^) = -ma{k,p), 



Xip, k- E) = Zip, k-E)+ j dq q' Zip, q; E) k; E) 
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with Z{p, q; E) given exactly by Eq. above. Note, in particular, that the homogeneous 
equation corresponding to (fTTj) requires no manipulation to be equivalent to that correspond- 
ing to Eq. ()B1|) . The relationship of X to the phase shifts can be deduced from Eq. ()B4|) 
and (jnni). It is: 



Re 



X{k,k;E) 



87 

k cot 6 

3m 



in agreement with Eq. 

In the case of nd scattering in the doublet channel we begin with the coupled equations of 
Ref. 0, 

which, again in the absence of a three-body force term, may be written in matrix 

form as: 

2 ? 

t{p,k)=Y{p,k) + - dqq^V{p,q;k)B{q;k)t{q-k), (B9) 

71 J 



with: 



and 



ttt tfs 
tst tss 



D 



Dt 




V = — In 

Apq 

Here, to leading order in 'jR, 

Dt{q; k) 

Ds{q\k) 



p^ + q^ + pq — mE 



p2 _^ q2 _ pg _ ^jt; 













;) 




^3 





-7s + xh' + Uf-k') 



(BIO) 

(Bll) 

(B12) 
(B13) 



To obtain Eqs. is now very simple. We just define: 

X = -m( ] (B14) 

Using Eq. ()B8|) to determine the relationship to the doublet phase shifts and K as in Eq. (fT^ 
we find that the t-matrix of Bedaque et al. should obey: 

3 1 



"'^ ' 8'j k cot S — ik' 
in agreement with Eqs. (12) and (13) of Ref. 3^. 
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equation for a rank-one separable two-body potential. 
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